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ABSTRACT 


Several techniques for calculating dose rates in the study of shelter entranceways 
are reported herein. Calculations have been verified by experiments. These calcula- 
tions alsu determine neutron energy spectra, a parameter that is difficult to measure 
and which is very useful in gaining insight into the design of protective structures. 
Roth the neutron flux spectra and dose rates are calculated for various positions along 
the centerline of a two-legged duct. Calculations were performed utilizing the 
ADONIS computer code which solves the Neumann approximation to the transport 
equation by Monte Carlo techniques. These Monte Carlo calculations were compared 
to second-leg calculations utilizing a semiempirical albedo formula, and were also 
compared to experimental measurements made with a tissue equivalent dosimeter. 
Included is a discussion of UNIGEOM, a generalized computer program for ADONIS 
qeometry, designed to be run on the IBM 7094, 
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teaults obteined by those who heave applied the information. 
This work sponsored by the Defense Atomic Support Agency. 


Pol ; 
4: es ee De ‘a 
: ‘ ‘ a 

} a al ie & ade : 

t as Pp) ae. ks ? et of is 

ry : Seat * a _ * 

£z 

Bor 


i 


nee eines Sheraton 


BepGidere smu a 


CLEARINGHOUSE FOR FEDERAL SCIENTIFIC AND TECHNICAL INFORMATION, CFSTI 
INPUT SECTION 410.11 


LIMITATIONS IN REPRODUCTION QUALITY OF TECHNICAL ABSTRACT BULLETIN 
DOCUMENTS, DEFENSE DOCUMENTATION CENTER (DDC) 


Ab 6/657) 
Cl |. — AVAILABLE ONLY FOR REFERENCE USE AT DDC FIELD SERVICES. 
COPY IS NOT AVAILABLE FOR PUBLIC SALE. 


By 2. AVAILABLE COPY MLL NOT PERMIT FULLY LEGIBLE REPRODUCTION. 
REPRO OUCTION WILL BE MADE IF REQUESTED BY USERS OF DOC. 


Bs A. COPY IS AVAILABLE FOR PUBLIC SALE. 
CI B. COPY IS NOT AVAILABLE FOR PUBLIC SALE. 
Cj 


3. LIMITED NUMBER OF COPIES CONTAINING COLOR OTHER THAN BLACK 
AND WHITE ARE AVAILABLE UNTIL STOCK IS EXHAUSTED. REPRODUCTIONS 
WILL BE MADE IN BLACK AND WHITE CNLY. 


Cr, ae 
TSL-121-2/65 DATE PROCESSED: 2 jem ieee 
PROCESSOR: ~ f). Nara. Le 
‘3 cr 


Cm! 


INTRODUCTION 


The U. S. Naval Civil Engineering Laboratory has for some time experimented 
with the streaming of nuclear radiation through the ducted entranceways of scale- 
model shelters, using 14.5-Mev neutrons as source particles. This energy approximates 
that of the prompt neutrons released in a fusion-type detonation, 

There are many techniques for calculating the dose rate of streaming neutrons, 
and this report compares results of the ADONIS Monte Carlo calculations with albedo 
calculations. All of the calculations shown have been verified by experiment. The 
great advantage in employing calculations is that costly and time-consuming experi- 
ments are minimized. 

The work reported herein is for the case of neutrons streaming through 3x 3-foot- 
square and 11x 11-inch-squore two-legged concrete ducts. The 14.5-Mev neutron 
source was placed just inside the duct entrance. 


ADONIS COMPUTER PROGRAM 


ADONIS! is a Monte Carlo computer program for solving the Neumann 
approximation to the transport equation governing neutron and gamma-ray streaming. 
It is written in FORTRAN II, version 2, for the IBM 7090, and ten tapes are required. 
ADONIS treats a problem having up to a maximum of 80 regions which, when com- 
bined, must cover all XYZ space where each region is a parallelepiped. It is 
important to note that the geometry diagram | in Appendix C of the ADONIS manual! 
contains an error in the designation of tl.e coordinate system. Also, it is necessary 
that the exit from the duct be in the positive Y direction and into the highest num- 


bered region. 
The particle flux is divided into energy groups called bins, of which a maximum 


. of 10 is permitted. For a neutron problem, appreximately seven nuclides may be used; 


however, for a gamma-ray problem more nuclides are possible, These nuclides must 
completely describe the material composition of each pcrallelepiped. The composition 
of each region may differ, but the composition within a particular region is limited to 
a combination of, at most, five elements. 

The geometric input required for ADONIS is time-consuming and complicated 
to prepare, even with the aid of various computer programs designed to simplify this 
task. In general, the geometry must be recomputed for each change in duct configu- 
ration. In many cases all that wiil be changed are dimensions such as length and 
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width of the legs, height of the duct, and the amount of material surrounding the 
duct. In order to facilitate the use of ADONIS in investigating the effects of such 
changes in duct parameters, a computer program was written which will accept the 
physical characteristics of the duct and as output will yield punched cards ready for 
an ADONIS run, These cards consist of the bounding plane cards (BPO), adjacent 
region cards (IRX}, complex plane cards (CMX), and horizontal and vertical plane 
cards (H and V). They define the geometry of ADONIS, The 8PO cards give the 
coordinates of the maximum and minimum extent of the regions along the axes, The 
IRX cards define which of the regions are back to back. The CMX cards define 
which of the regions back to more tiian one region. 

Although subroutines called EASYGEOM and TEST G are available to 
facilitate preparing the input cards to ADONIS, mistakes are easily made in geome- 
tries as complex as those associated with duct problems. To remedy this situation, a 
computer program, UNIGEOM, was written at NCEL. The source deck of this program 
is listed in Table I. This program incorporates a modified EASYGEOM program as 
the main routine and TEST G end its subroutines as a subroutine. The bounding plane 
cards required as input to EASYGEOM are calculated from the physical parameters 
of the duct by means of a subroutine: The UNIGEOM program all. ws three such 
subroutines. The subroutines included in this report are for a two-'sgged duct and 
three-legged coplanar and non-coplanar ducts, Other configurations may be incor- 
porated by substituting "duct" subroutines and appropriately changing the calling 
" sequence in the main program. 

ADONIS considers the following neutron interactions: isotropic and anisotropic 
elastic scattering, discrete and continuous inelastic scattering, scattering by hydrogen, 
and the (n, 2n) reaction in beryllium absorption. The gamma portion of the ADONIS 
program considers the Compton scattering of photons where the photoelectric effect 
and pair production are regarded as absorption. Particle histories in ADONIS are 
terminated either when an absorption occurs, the particle is degraded below an 
arbitrary cut-off energy, o kill occurs due to splitting, or the particle escapes. Cross- 
section data used in ADONIS are called from an organized-element-Jata tape. A 
source-data tape containing the position coordinates, direction cosines, and the 
initial energy of the neutron or gamma-ray source is required in addition to the 
organized-element-data tape and the geometry-input data. 

The source-data tare is generated for the case of either the point source 
incident radiation of neutrons or gamma rays by PSGEN, a rrogram modified to be 
anisotropic; in the case of gamma rays produced by neutron interaction, the source- 
data tape is generated by the GASP! program. The input to PSGEN, which is modified 
to be a 2 source, describes the source as the probability that part'cles exist below a 
particular energy. The energy is evaluated for probabilities in 0.1 increments of 
probability from 0 to 1.0. ADONIS is called in through the IBM 7090 monitor system. 
ADONIS contains dump and edit procedures, selected by turning sense switch 1 "on" 
at the time these procedures are to be initicted. After dumping, ADONIS may be 
restarted. Turning sense switch 2 "on" prints on-line the particle history numbers; 
turning sense switches 3 and 4 "on" writes an off-line trace of each particle path on 
tape A3. 


Particle splitting is employed to improve the efficiency of the calculation by 
assigning importance weights to each of the regions. Regardless of the importance 
selected for each region, if the ADONIS program is allowed to run for o sufficiently 
long time, the answers furnished by it will be approximately the same. However, in 
making many calculations it has been determined that the proper selection of weighting 
functions can significantly decrease the amount of running time required. The optimum 
value af imnortance for each region may best be selected by a semiempirical approach. 
First, the geometry is studied and a weight, inversely proportional to imporicnce, is 
assigned to each region based upon the best knowledge of neutron or gamma-ray 
streaming available to the investigator. Using these values, the program is run for 
a short period of time, after which the calculated flux in each region is reviewed, 

The weight of the region is then adjusted in a way which will enhance the probability 
of neutrons going in the desired direction while reducing the probability of neutrons 
going in undesired directions. 

It has been experienced that approximately five tries are necessary to optimize 
the weight for each region in a duct with two legs and approximately 10 tries for 
similar optimization in a duct with three non-coplanar legs. It shouid pe emphasized 
that a proper selection of weight allows a greater percentage of computer time to be 
devoted to following particles of interest; thus, running time for statistically signifi- 
cant determinations of flux may be considerably reduced. A proper selection of 
inportance functions and an adequcte running time may be judged by comparing the 
standard deviation of particle flux to the flux for each energy bin and dust iegion. 
ADONIS calculates the flux and particle dose in each region. The atter is accom- 
plished through either o fine or coarse response function. These two functions differ 
only in that for the fine function, the response is evaluated at the same energies used 
for cross sections on the element-data tape, and for the coarse function, the response 
is averaged over ecch energy bin. Also, the standard deviation of dose is evaluated 
only where fine response functions are used. The response function is simply a function 
that converts particle flux into dose such as relative biological effectiveness (RBE) 
for the case of neutrons. 

ADONIS also provides, as output data for each region, the number of absorptions, 
' births, degradarions, and deaths; flux per eneray bin; deviation of flux per energy bin; 
total number of reflected and transmitted particles; particle balance; absorption by 
element and region; and coarse dose per energy region. The neutron fl.x is reported 
as the track length per source particle per unit energy in energy interval j per unit 
volume for finite regions. In the case of regions infinitely bounded, the total flux is 


given. 


Analysis of ADONIS Results — 


When ADONIS is restarted, each run prints only the flux and dose associated 
with that run, rather than the best estimate based upon total running time. It was, 
therefore, necessary to calculate these parameters as follows: 
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RBE.(E) 
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best estimate of dose based upon n runs 
energy increment of jth bin 

number of histories of it? run 

number of runs 

number of energy bins 


relative Lological effectiveness averaged over the energy 
range of the j'" bin 
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standard deviation of dose for n! run 
best estimate of standard deviation of dose based upon n runs 


standard deviation of flux for i bia 


best estimate of flux based upon n runs 


where D,, is the incident dose due to the total number of particles emitted by the 
analytically simulated neutron source, Since the analytical treatment considers 

these to be emitted in one hemisphere, the value of incident dose is divided by 27, 
the number of steridiacs ina hemisphere, In actual practice, ADONIS calculates 

the dese per incident source particle, and thus the associated value of incident dose 

is one source porticle multiplied by the RBE of the source particle, This is numerically 
equal to the RBE of the source particle energy. 


DESCRIPTION OF EXPERIMENTS 


In experiments at NCEL, 14.5 Mev neutrons approximated the energy of neutrons 
released in a fusion-type detonation. Neutrons of this energy were available from a 
neutron generator utilizing the T(d,n)He4 reaction. A two-legged duct ° x3 feet and 
a three-legged duct 11x11 inches were constructed for the experimental phases of 
this work, The latter was constructed in such a way that the orientation of the third 
leg could be altered to form either the coplanar or non-coplanar case. Plan and 
elevation drawings of these ducts, including dimensions, will be found in Figures 1, 

2, and 3, The neutron source was approximated as a point source for the purpose of 


caleulations. In actuality, the active area of the target was a circle approximately 
7/8 inch in diameter, Although this target radiated more or less isotropically, atten- 
tion was taken to orient the geometry in such a way that the plane of the target was 
located concentric to the centerline of symmetry and perpendicular to the duct 
approximately one centimeter from the end of the first leg and inside the duct. The 
same location was used for theoretical calculations, These calculc tions were made 
for both the 3-foot and the 11-inch duct. However, experimental neasurements 

- were made only with the 3-foot duct, since the physical size of the dosimeter 
significantly perturbed the flux and in certain cases simply would not fit. 

There is some present concern for the statistical variation of the experimentally 
measured dose rates, the calibration of the instrumentation utilized, and the determi- 
nation of tne inciden? dose, D,. Nonetheless, it is believed that these results are 
certainly accurate to within an order of magnitude. Although pluns exist for 
repeating these measurements with more care and augmenting them with measurements 
performed with helium and argon-methane detectors and foil activation techniques, 
this will not be done for some time. Such work really has little bearing upon the 
comparison of ADONIS calculations reported herein, since it will, in all likelihocd, 
simply add more significance to the determinations. 

While reviewing the results of Doty's? experiments, Songs observed thot a 
semiempirical albedo formula could be utilized to fit the duta. He utilizes Monte 
Carlo results reported by Ballistic Research Laboratory4 in order to determine the 
constants of his formula. He then calculated the dose in the second leg of the 
3-foot duct and found it to be in good agreement with Doty's experimentally 
~ measured values. While it is que that many different functions could undoubtedly 
fit Doty’s data, a physical explanation of neutron streaming can be easily made 
utilizing albedo concepts. This was the motivating reason behind the albedo devel- 
opment for neutron streaming through ducts. For additional volidity, Song calculated 
the dose in the second leg of the il-inch duct, and the ADONIS program was 
utilized to determine the dose in both ducts by Monte Carlo calculutions of the 
Neumann approximation to the transport equation, 


MONTE CARLO CALCULATIONS 


"The program UNIGEOM (Table !), was written to calculate ADONIS input 
geometry from the physical parameters of the three duct configuratons shown in 
Figures 1, 2, and 3. Since TEST G is incorporated in this program, the physical 
dimensions of the duct are checked for consistency, and the same diagnostics are 
written on the output tape as would be had EASYGEOM and TEST G been run as 
they were originally intended and not incorporated in the present program, 
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Figure 1. Two-legged duct with regions and physical parameters. 
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Figure 2. Three-legged coplanar duct with regions and physical parameters, 
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The appropriate geometry subroutine is indicated by a value on the first input 
card, Input data are read by the main program from tape A2 and are communicated 
to the proper routine through Common, Input data in the format F 10.5 are read in 
sequence XL], XL2, XL3, WI, W2, W3, TSI1, TS12, TS21, TS22, TS31, 1532, TT, DI, 
D2, 03, DA, DB, DC, DAI, DA2, DA3, TA, TB, TC, TD, HH, H1, H2, XX, Table It 
lists the values of these parameters for each of the three duct configurations. In 
some cases these values are zero, but for proper sequencing must be provided, All 
dimensions are given in centimeters to be consistent with ADONIS. The various 
regions in the duct are also identified in Figures 1, 2, and 3, It will be noted that 
a matrix of nine regions cuts each important leg in the duct. The center region of 
this matrix is an area } centimeter square and may be utilized as a detector region 
useful for comparing calculations to experimental dosimetry ata point. The geometry 
coordinates of the subroutine calculations are written as X(A, B), where A has the 
value given in the following table. ) 


Coordinate Coordinate Number 
-X 1 
+X 2 
-Y 3 
+Y 4 
a7 5 
+Z 6 


Coordinates are indicated by X(A, B), where A is the 
coordinate indicator ard B is the region number, 


This table also compares with values in ADONIS notation. The region number 
is given by B, The output of UNIGEOM is written on tapes A3 and B5, the latter of 
which is punched as geometry cards required for ADONIS input. First the BPO cards 
are punched, as identified in the preceding table, one card for each region. These 
coordinates are punched in E 14.5; next, the cards for the coordinates of each com- 
plex plane for both horizontal (H cards) and vertical (V cards) are punched. A 
complex bounding plane is created for the surface of one region if it is adjacent to 
the surfaces of more than one other region, One or more cards is devoted to the 
coordinates of the horizontal and vertical fields of each complex bounding plane. 
The adjacent region (IRX) cards and the complex plane (CMX) cards needed for 
ADONIS input are also provided by the program. Should there be an error in 
geometry or physical dimension that prevents the calculation of geometry, various 
diagnostics are provided. ADONIS input cards calculated by UNIGEOM for the 
parameters of Table II are listed in Table II}. This information is provided as 
sample problems. 
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COMPARISON OF THEORY TO EXPERIMENT 


Figure 4 shows the neutron energy flux calculated by ADONIS in several 
regions associated with the 11-inch two-legged duci, normalized by D,/21. These 
regions may be identified by referring to Figure 1. In Figure 4 it is seen that the 
flux, due to source neutrons, steadily decreases with distance along the first leg. 

A few source neutrons are also found in region 14, the region in the second leg 
adjacent to the corner. These neutrons must have arrived by direct penetration, 
grazing the inside comer of the duct. Had it been otherwise, their energy would 
have been degraded because of interaction. No source neutrons are found at other 
regions in the second leg. 

It should be noted that the flux approximates a straight line on log-log graph 
paper for energies lower than|~1 Mev. In order to produce the results shown in 
Figure 4, it was necessary to pperate the 7090 computer for approximately 23 hours. 
Because of this long operating time and the straight-line relationship noted, it was 
decided that in the future ADONIS. runs the energy of the low-energy cut-off point 
would be raised. This should resuli in approximately the same statistical deviation 
with a much shorter running time. 

The dose, corresponding to the flux shown in Figure 4, is shown in Figure 5 
for the same energy regions. Here, the dose has been normalized as previously 
described, by dividing with D,/21. It should be noted that the dose in the first leg 
falis off as an inverse square function of distance, and in the second leg the dose 
falls off at a much steeper rate. This corresponds to similar results obtained during 
gamma streaming experiments. In Figure 5 are shown the values of dose calculated 
utilizing Song's semiempiricai albedo formula. It is seen that the agreement is good. 

Also shown in Figure 5 are the Monte Carlo/Albedo calculations at two lower 
source-particle energies, utilizing the Oak Ridge National Laboratory code. It is 
seen that the ORNL results do not agree too well and are consistently high. It should 
be noted that the assumptions of corner lip secondary source preclude utilizing the 
albedo formula for the corner region and makes the results of the region adjacent to 
the corner somewhat questionable. In Figure 6, the value of D/D, has been corrected 
to remove the effects of distance. This has been done by multiplying the normalized 
dose values shown in Figure 5 by the factor 4” where 4 is the length .i the duct leg, 
and the exponent, n, is determined from the linear slope of the dose tines shown in 
Figure 5. If the duct were ideal and long, the figure resulting from this new calcula- 
tion would ideally be a step discontinuity occurring at the comer region. Because 
of the small length for the duc! (~140 centimeters), this idealized shape is not realized 
and the resultant S-shape is obtained. 

In the preceding figures are seen the first collision, multiple collision, and rad 
tissue equivalence dose rates. In Figure 7, the flux spectrum resulting from ADONIS 
calculations is shown for the 3-foot duct where the leg length is longer than that for 
the Il-inch duct. Because of incr used computer running time required for the larger 
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duct, the low-energy cut-off limit was increased to | kev and the computer running 
time was considerably decreased. This resulted in an increase in the stordard devia- 
tion of the flux, as shown by the bounds of each of the bars comprising the histogram. 
This flux is very similar to that shown for the small 11-inch duct, and the slopes of 
the spectrum lines below 100 kev approximate those shown in Figure 4 between 

100 kev and 1] Mev. There is considerable difference between the spectra shown for 
the large duct in Figure 7 and those shown for the small duct in Figure 4. It is felt 
thet this could be explained by the fact that the increased length of the large duct 
affords a greater opportunity for neutron interaction and, thereby, degradation of 
energy spectra. Results from the two ducts are, however, similar. 

In Figure 8 is shown the dose of the !arge duct, normalized as before, the dose 
calculations utilizing the semiempirical albedo formula of Song, the ORNL dose cal~ 
culations, and the experimentally measured dose obtained by Doty, It is seen that 
Song's and Doty's work are in good agreement; yet their dose measurements are 
greater than those calculated by the ADONIS code and less than those calculated 
by the O5R ORNL code. It should be noted that up to this point the importance 
functions of the regions of the two-legged-duct problem were selected to eliminate 
most corner-lip diffusion. 

Another Monte Carlo ADONIS calculation was performed, keeping all 
conditions similar to those in the previous calculation of the large duct, with the 
exception that the importance of the rectangular region between the two legs was 
now selected to allow diffusion over the entire rectangular region. The results of 
these calculations are shown in Figure 8 by the dashed lines. It is seen that these 
results lie somewhat above the results of Song and Doty and agree very well with 
ORNL calculations. 

In Figure 9 is shown the dose corrected for fall-off because of length. This is 
done by multiplying the dose of Figure 9 by 4" as before, Here, the effect of the 
comer is more prominent, and the resulting curves more closely approximate the ideal 
step function. It thus appears that it is necessary to carefully consider the size of the 
comer lip secondary source for albedo calculations and, indeed, possibly include some 
diffusion effects. A great deal of caution shouid be exercised in fitting formulas, 
because when tne standard deviations of the data are large —as found both experi- 


. mentally and theoretically —it is usually possible to obtain what appears to be a 


good least-squares fit with almost any curve. On the other hand, it also appeors that 
it is incorrect to assume maximum diffusion through the rectangular region; rather, it 
is believed that a more involved type of interaction is occurring between the two 
effects and that such interaction could be approximated by selecting the appropriate 
importance function for the rectangular region adjacent to both legs. Again, calcu- 
lations which would erily this conclusion were not within the scope of the present 


work, 


12 


genset es Sill Sh 


ba Beckie GRAS gl we hee OE: A ete 


~~ : 
_ 
7 
~~ 


on ow a oe em one De 


Region Numbers 
2 (er 6 Cece 


15 


14 
24 


“tar l6 
1 
: 25 


10°!7 


10 102 102 104 
Energy (e 


Figure 4, 11 x 11-inch duct 
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Figure 5. 11 x 11-inch duct, normalized neutron dose. 
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CONCLUSION AND RECOMMENDATIONS 


The Monte Carlo calculations of neutron streaming through model ducted 
entranceways, performed utilizing the ADONIS code, and the Monte Carlo/Albedo 
calculations utilizing the ORNL codé are in good agreement. Protection factors 
were found to be proportional to the duct cross-sectional area. It was found that 
albedo calculations of dose must be based upon a much grecter corner penetration 
than previously assumed. This could decrease previously reported neutron protection 
factors by approximately a factor of 10, It remains, however, for these conclusions 
to be verified by carefully controlled experiments. By performing analytical calcu- 
lations, as described, protection factors of structures may be readily determined with 
a fraction of the effort required by experimental methods, It is, however, necessary 
that all analytical techniques be verified by careful experiments. Once verified, 
the techniques may be freely applied in general extensions to the experiments. 

It is recommended that the physics of neutron penetration between the first 
and second legs by methods other than streaming through the duct be more intensively 
investigated in order to establish the validity of albedo representations. It is also 
recommended that the ADONIS and ORNL codes be accepted as preferred analytical 
techniques for determining neutron streaming through ducted entranceways. The 
ADONIS code is applicable to most problems. The ORNL code, while requiring less 
running time, is applicable only to the case of a duct through o dense, infinitely 
bounded medio, and thus it has limited practical applications. 
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Table Il. Physical Parameters of Ducts 


Subroutine 
Purametess L3Z 
(Two-Legged) 

Fisst Card 

Colunms 9 & 10 43 47 

Column 20% 1 2 
XL} 200 cma 167 em 
XL2 190 151 
XL3 0 14) 
wi 92 72.39 
W2 100 80.01 
W3 0 72.39 
TSH 29 16.51 
1812 35.6 27 .94 
7821 35.6 16.5] 
1822 35.6 35.56 
TS31 0 27 .94 
1§32 0 16.5] 
TT ! 0 
DI 41 10 
D2 41.0 9 
D3 0 0 
DA 1.0 0 
DB 1.0 0 
Ac 1.0 0) 
DAI 0 1.0 
DA2 0 1.0 
DA3 0 0 
TA 6) 61 
7B él 61 
Tc 6} él 
TD 6} 6] 
HH 60.9 60.96 
H1 1o.5 16.51 
H2 16.5 16.51 
XxX Q | 0 


Y Total number of regions 
2 Subroutine number 


J Duct dimensions 
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(Three~Legged Coplanar) | (Thiee-Legged Non-Coplanar) 
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